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Shocks and dust survival in nearby active galaxies: 
implications for the alignment effect 



o 
o 

(N 

< 
(N 



> 
O 

m 
oo 
o 



M. Villar-Martm 1 , D. De Young 2 , A. Alonso-Herrero 1 ' 5 , M. Allen 3 , L. Binette 4 

1 Dept. of Natural Sciences, University of Hertfordshire, College Lane, Hatfield- Herts, ALIO 9AB 

2 National Optical Astronomy Observatory, 950 North Cherry Avenue, P.O. Box 26732, Tucson, AZ 85726, USA 

3 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, Maryland 21218, USA 

4 Institute de Astronomia, UNAM, Apartado Postal 70-264, D.F. 04510, MEXICO 

5 Present address: Steward Observatory, The University of Arizona, Tucson, AZ 85721, USA 



ABSTRACT 



45 

6 
to 

CO 

> \ 1 INTRODUCTION 

> I ■ The discovery of high levels of polarization in the extended 
C\3 continuum of many powerful radio galaxies at redshift z >1 

(Cimatti et al. 1997; Vernet et al. 2001) proves that the 
alignment between the optical and the radio structures (Mc- 
Carthy et al. 1987; Chambers et al. 1987) observed in most 
powerful radio galaxies has an important contribution of 
light from a hidden quasar scattered towards the observer by 
dust and/or electrons. Several reasons support that dust is 
the scattering agent, rather than electrons (see Cimatti et al. 
(1998) and De Young (1998) for a discussion). This interpre- 
tation faces some problems: there is now increasing evidence 
for the existence of fast shocks in the extended gas of distant 
radio galaxies generated during the interaction between the 
radio jet and the ambient gas (e.g. Best, Rottgering, Lehnert 
2000). De Young (1998) raised the question of the unlikely 
survival of dust under these extreme conditions. The author 
calculated the survival of dust grains for a wide range of 
parameters appropriate for the conditions expected in dis- 
tant radiogalaxies. De Young found that for most of the 
cloud configurations the grains are destroyed as a scattering 
population by sputtering processes in a time much less that 
minimum radio source lifetime of ~10 million yr. 



One of the most popular explanations for the so-called alignment effect in high 
redshift (z >0.7) radio galaxies is the scattering by dust of the hidden quasar light. As 
shown by De Young (1998) a problem with the dust scattering model is that the short 
destruction time-scale for dust grains means that they will not survive the passage of 
the radio jet. 

We investigate the survival of dust in the extended ionised gas of nearby active 
galaxies with jet/gas interactions. We discuss the implications on the alignment effect 
of high redshift (>0.7) radio galaxies. We conclude that although shocks are likely to 
destroy dust grains in regions of interaction, dust might survive in enough quantities 
to scatter light from the active nucleus and produce alignment between scattered light 
and the radio structures. 

We propose an observational test to investigate the existence of dust in shocked 
regions based on the sensitivity of calcium to depletion onto dust grains. 
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We investigate in this paper the existence of dust in 
the extended gas of nearby active galaxies where jet/cloud 
interactions are occurring. There are several low redshift ac- 
tive galaxies where evidence for shocks (originated in jet/gas 
interactions) and dust in the extended gas has been found. 
Cygnus A is an example: dust illuminated by the active nu- 
cleus is observed in optical linear polarization maps which 
show a giant bipolar reflection nebula (Tadhunter, Scarrot & 
Rolph 1990). Tadhunter (1991) suggested that the compo- 
nents of high velocity gas discovered in the extended narrow 
line region could be the result of the interaction between 
the radio jet and the gas. Another example is PKS1932- 
46. Villar-Martm et al. (1998) proposed that the extended 
gas is interacting with the radio structures. The existence 
of dust in the extended gas is suggested by the lack of 
[Call] AA7291, 7324 emission (see Fig. 4 in Villar-Martm et 
al. (1998)). The absence of these lines is likely to be due to 
the depletion of calcium onto dust grains (§5). 

Therefore, jet/cloud interactions and dust exist in the 
extended gas of some radio galaxies. However, in these spe- 
cific cases it is not clear whether dust is mixed with the 
shocked gas. Grain destruction might have occurred in the 
shocked gas, but grains could still exist in the non perturbed 
regions. 
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We discuss in this paper the case of two active galaxies 
where the apparent coincidence between the radio and dust 
rich structures strongly suggests that the shocked gas is 
mixed with dust: the Seyfert 2 galaxy NGC1068 (Seyfert 
1943) (§2) and the radio galaxy PKS2152-69 (Westerlund & 
Smith 1966) (§3). By applying models appropriate for the 
shock conditions in NGC1068, we investigate whether dust 
grains are likely to survive the passage of the radio jet. In 
§4 we propose an observational method to search for dust 
in the shocked (or non shocked) gas of low redshift active 
galaxies. Discussion and conclusions are presented in §5. 



2 NGC1068 

2.1 Synchrotron emission from the radio jet? 

The Seyfert 2 galaxy NGC1068 (Seyfert 1943) has been 
widely studied at all wavelengths (see for instance Astro- 
physics and Space Science, v. 248). Mid Infrared (MIR) ima- 
ges reveal a remarkable similarity between the MIR and the 
radio structures, suggesting that they are spatially associa- 
ted (Alloin et al. 2000, AL00 thereafter). If this is the case 
and if the MIR emission is produced by dust, then NGC1068 
could be an example of an object where dust has survived 
the passage of the radio jet. 

The authors identified four structures (in addition to 
the core) in NGC1068^ at 11.2 and 20.5 jum, labelled as a, 
b, c and d (we adopt the same nomenclature). Regions a 
(1.8" X4.0") and b (1.3" x 1.3") are located along PA=35 deg 
to the northeast at mean distances from the core of 3.7" and 
1.8" respectively. Regions c (1.3" xl. 6") and d (1.3" xl. 3") 
are along PA=210 deg to the southwest at mean distances 
from the core of 2.1 and 3.5 " respectively. 

First we have investigated the possibility that the MIR 
emission is due to synchrotron radiation from the radio jet, 
rather than thermal emission from dust. This is suggested by 
the close correspondence between the radio and MIR struc- 
tures. Jets are rarely observed at MIR and optical wave- 
lengths but some cases exist, such as the radio galaxy M87 
(e.g. Stocke, Rieke & Lebofsky 1981). For NGC1068 we have 
calculated the expected synchrotron emission at 20.5 /zm 
(the emission will be fainter at 11.2 (im) from the radio jet 
in two of the bright MIR structures as identified by AL00 (b 
and c). We measured the surface brightness at 4.9 GHz in b 
and c using the radio map of Wilson and Ulvestad (1987). 
Assuming S oc u~ a and using the spectral index between 
4.9 and 15 GHz measured by these authors in both regions 
b (a=2.1) and c (q=1.2), we calculated the expected syn- 
chrotron emission at 20.5 /jm. These values are upper li- 
mits, since the spectral index of the synchrotron spectrum 
is often observed to steepen with frequency (e.g., Cyg A. 
Kellermann, Pauliny-Toth & Williams 1969). The results 
(presented in Table 1) show that synchrotron emission does 
not make a significant contribution to the MIR emission. 
We therefore conclude that the MIR extended emission in 
NGC1068 is due to dust (as claimed by e.g. AL00). 



* Throughout this paper we will adopt a distance to NGC1068 
of 14.4 Mpc where 1" corresponds to ~72 pc 



2.2 The survival of dust grains 

Have the dust grains survived the passage of the radio jet 
in NGC1068 ? Is a mechanism of replenishment required to 
explain the existence of dust mixed with the shocked gas? 
De Young (1998) studied the response of dust grains to the 
passage of the strong shocks associated with the radio source 
in an environment where the dust is located in gas clouds of 
high density contrast and small filling factors embedded in a 
more diffuse medium. He explored a wide range of parame- 
ters. The factors that influence the evolution of the grains 
are the initial cloud gas density, the cloud size, the initial 
grain size, and the speed of the shock associated with the 
radio jet. Indirect factors arise from the density of the inter- 
cloud medium and the presence or absence of magnetic fields 
in the clouds. We have run similar models with parameters 
appropriate for NGC1068 and we have investigated whether 
dust can survive the passage of the radio jet in this galaxy. 



2.2.1 Model input parameters 

The parameters for NGC1068 are the following: 

• Initial cloud gas density (n w ): Spectroscopy of the re- 
gions not interacting with the radio jet implies densities of 
100 cm' 3 [Axon et al. 1998]. Since this gas is not affected 
by the shocks, we have considered this value as the initial 
cloud density, before the passage of the radio jet (the mea- 
sured density of the clouds in the region where the jet is 
interacting with the gas is 10 4 ' 5 [Axon et al. 1998]; in our 
interpretation, such high densities arise naturally from shock 
compression) 

• Number density contrast nh/n w : (where nj, is the den- 
sity of the intercloud medium). The measured temperature 
in the ionised gas is T e ~1.15x 10 4 K (Axon et al. 1998). We 
expect the hot phase temperature to range from 10 6 to 10 7 
K (e.g. Forman, Jones & Tucker 1985, Shirey et al. 2000, 
Netzer & Turner, 1997). If the two phases are in pressure 
equilibrium, this implies that the number density contrast 
fin /n w (which is the only intercloud medium parameter that 
matters for grain evolution) is in the range 10" 2 to 10" 3 . 

• Speed of the shock associated with the radio jet v s : 
Axon et al. (1998) report a velocity shift between the split 
components in the optical emission lines of 1500 km s" 1 at 
the position of the interaction between the radio jet and the 
ionised gas. We assume that the velocity of the clouds accel- 
erated by the shocks is ~750 km s _1 (Axon et al. (1998) de- 
duce shock velocities of at least 700 km s _1 ). In the absence 
of entrainment processes, the cores of the warm clouds will 
be accelerated to v w — v s (ph/ pw) ' 5 ~ v B (nh/n w ) ' 5 , where 
v w is the final velocity of the warm clouds and v s is the 
shock speed in the hot phase. This implies shock velocities 
in the hot phase in the range 7xl0 3 km s _1 and 2xl0 4 km 
s _1 , depending on the number density contrast. 

• Clouds size (Rdoud). as the uncertainty on this value 
is high we have considered as upper and lower limits for 
the radius of a cloud as 1000 pc and 1 pc respectively. We 
find that the only time the initial cloud radius makes any 
difference is in the B=0 (no magnetic field) case, which is 
rather unlikely (see below). 

• Dust properties: we have assumed (as AL00) the MNR 
distribution n(a) = K a a~ 3 5 , where n{a) is the number den- 
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Or> 20.5 
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1.8 
2.1 


1.3x1.3 
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0.031 
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2.1 
1.2 


4.3e-7 
6.8e-8 


0.65 
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Table 1. Comparison between the predicted surface brightness at 20.5 /an due to synchrotron emission (SBjq^' 1 ) and the observed 
surface brightness (SBJ^y"™ 3 ) at 20.5 lira (ALL00). D is the distance to the core in arcsec. Ext is the projected size. SB4 .qgHz is the 
surface brightness in b and c at 4.9 GHz. a is the spectral index measured by Wilson and Ulvestad (1987) between 4.9 and 15 GHz. 
Synchrotron emission makes a negligible contribution to the MIR in NGC1068. 



sity of grains and a is the grain radius. Upper and lower 
limits on a were taken as 1000 and 10000 A respectively . 

2.2.2 Models results 

We show in Table 2 the results of the models for the range 
of parameters appropriate for NGC1068. The final column 
indicates the survival time of dust grains of size i?o(A) under 
the conditions considered in that specific model. (Grains are 
considered to be destroyed when the grain radius falls below 
one percent of the initial radius, leaving only 10 -6 of the 
initial mass). For the models in the last three rows, the scale 
for dust destruction is >10 6 yr. 

Therefore, grains of 1000 A radius are destroyed in less 
than 1000 yr years for the shock speed and density contrast 
values expected in NGC1068. Only if there is no magnetic 
field in the cloud, or in the case of 'bare grains' (see below) 
can small grains survive more than a million years. The large 
grains (10 4 A ) can survive in low velocity shock conditions 
for very long times. Since almost all the grains are at the 
small end of the scale in the MNR distribution, ~1000 years 
after the passage of the radio jet most dust should have been 
destroyed. 

If the outer part of the radio jet interacted with the 
gas in region 6, when did this interaction happen? Using 
an advance jet speed in the range O.lc-O.Olc (e.g., Scheuer 
1995)[] and the distance between the outer edge of the jet 
and region 6 ~4 arcsec (288 pc in projection), the interaction 
should have happened more than 10 4 yr ago. Therefore, ac- 
cording to the model results the dust should have been des- 
troyed. However, there is no evidence for dust destruction 
in NGC1068. 

We propose several explanations: 

• Dust has not been destroyed 

- The outer radio source did not interact with the gas 
in regions b and c. On the contrary, the interaction could 
have happened very recently (it might still be going on) 
with younger radio structures and there has not been 
enough time for dust destruction. 

An upper limit for the time elapsed since the interac- 
tion can be calculated from the cooling time of the 
shocked gas. If the shocked gas is still hot enough 

t Based on the similarities in radio source size and velocity be- 
tween Mrk 3 and NGC1068 Capetti et al. (1999) proposed an up- 
per limit for the age of the radio source in NGC1068 of ~1.5xl0 5 
yr. Using the projected distance between the outer edge of the ra- 
dio lobe and the core (~6 arcsec) the jet advance speed is >0.01e) 



to produce ionizing photons (as proposed by Axon et 
al. (1998), see also below), the time elapsed since the 
interaction should be shorter than the cooling time. 
For shock velocities in the range 500-1000 km s -1 , the 
shocked gas will cool to the recombination tempera- 
ture of hydrogen in a time-scale t ra d,& Myr given by 
trad,6 ~ 1.9 n' 1 v W 3 where v W 3 is the shock velocity 
in the clouds in units of 1000 km s -1 (Bicknell, Dopita 
& O'Dea 1997). For «„=700 km s -1 and n m =100 cm -3 , 
the shocked gas will cool in ~6.7xl0 3 yr. The interac- 
tion happened, therefore, less than 6.7xl0 3 yr ago. This 
might point to a recent interaction, but it is not possible 
to assure that it happened less than 1000 years ago. 

- The magnetic field in the cloud is B—0. In this 
case grains survive the slower shock speeds, since the 
cloud cools first and there is no betatron effect (Spitzcr 
1941). If B=0 and one turns on saturated conductivity, 
the grains still survive a long time. The cloud heats, 
and sputtering increases, but as the cloud heats it also 
expands, and the lower number density greatly reduces 
the sputtering rate. The small grains are still at 0.73 of 
their original radius after 2xl0 4 yr. The case of null mag- 
netic parameter is, however, unlikely. In all the galaxies 
we know about, the interstellar medium is permeated 
with a magnetic field. For active galaxies in particular, 
where the presence of synchrotron radiation implies a 
magnetic field, one would expect a significant magnetic 
field in clouds near the nucleus. We therefore reject this 
interpretation. 

- The clouds consist of an external shell of ionised 
gas and denser cores where the dust has been protected 
from the intense UV radiation from the active nucleus. 
After the radio jet has hit the clouds and the shocks 
have crossed the outer warm ionized layers (with n ~100 
cm' 3 ), the inner dusty cores become exposed. For the 
cores, rih /w^ ore could be so low that the shocks would be 
very slow and the grains could survive (n^° re is the gas 
density in the cores). This model has been proposed by 
Bremer, Fabian & Crawford (1997) to explain the align- 
ment effect in HzRG as due to scattered light from the 
hidden quasar, since an intense dust scattered continuum 
is expected along the radio axis. This is supported by 
the remarkable similarity between the radio structures 
and the J and H band polarized flux morphologies (due 
to dust scattered light) in spatial scales of ~5" from the 
nucleus (Packham et al. 1997) . The enhanced MIR emis- 
sion at the position of the radio structures also supports 
this interpretation (see §2.1.3). 

As an example, if the density contrast between the 
hot phase and the dense cores is nu /rC re = 10 -5 , a shock 
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Table 2. Grain survival rates for an internal cloud shock velocity of v w =750 km s 1 and the different parameters discussed in the text. 



propagating at vh~2x 10 4 km s _1 in the hot phase would 
advance with v!^ re —63 km s _1 within the dense cores. 
The destruction of the dust grains will be very inefficient 

with such slow shocks (De Young 1998). |^] 

For such density contrast, and assuming nh=0.1 cm -3 , 
the density of the dense cores is n^ re =10 4 cm -3 . If 
this phase is in equilibrium with the hot phase, then 
TS° re =10 K. These properties correspond to a molecu- 
lar phase. Lepp et al. (1985) showed that obscuration 
by gas and dust allows the existence of a stable (cooled 
by molecules) phase (T ~10 2 K) in the ISM heated and 
ionised by a soft X ray continuum. A column density of 
Air~10 22 cm -2 would be necessary to allow the exis- 
tence of this phase. If n„ = 10 4 cm -2 , this implies that 
the thickness of the 'protective' shell is ~0.3 pc, which 
is in good agreement with the minimum assumed total 
cloud sizes (1 pc). 

- In the case of bare grains (i.e. not embedded in 
protective clouds, but distributed uniformly trough out 
a homogeneous hot circumgalactic medium) the dust 
grains are very slowly eroded, with the smallest grain 
(1000 A) still being at 0.91 its original radius after 2xl0 4 
yr. A mechanism to generate this dust distribution is, 
however, unknown (De Young 1998). 

• Dust has been destroyed and a mechanism of replenish- 
ment is at work 

- Dust production in situ by jet induced star forma- 
tion is not viable, since there is no evidence for such 
a process in NGC1068. Another possibility is that dust 
grains are convected from the dusty interior of the galaxy 
outward with the continuously outflowing radio jet. The 
issue of grain survival in a hostile environment might, 
however, still be a problem (De Young 1998). 



In summary, the close association between the extended 
MIR and radio structures in NGC1068 suggest that the 
shocked gas is mixed with dust. Calculations of dust sur- 
vival for NGC1068 suggest that the dust should have been 
destroyed. There might be a mechanism that replenishes the 
dust grains, however, there is no evidence for it. Other possi- 
ble explanations for this inconsistency are a) the interaction 
happened recently and there has not been enough time for 



t Notice that the emission lines still require shock speeds of 
v w =750 km s — 1 and the shocks slow down to v^ re when they 
enter the molecular core. 



dust destruction b) the clouds contain inner molecular dusty 
cores where dust destruction is less likely to happen. 

2.3 The heating of the dust 

Why is the dust emission enhanced in the regions crossed 
by the radio jet? 

One possibility is that the AGN radiation field is highly 
anisotropic (Axon et al. 1998). An alternative explanation 
is that, in addition to the continuum from the active nu- 
cleus, there is another mechanism working along the radio 
axis that heats the dust locally to higher temperatures, en- 
hancing the radiation. 

The shocked gas can heat the dust locally via two diffe- 
rent mechanisms: impact with the dust grains and radiative 
heating. 

The most efficient heating by impact is due to electrons. 
In this case, the temperature of the dust grains is given 
(provided that T gaa » T dust ) by: 

_ 3 k T gas 3/2 U m e £ 1/4 
Jdust — [( — j 

m e I a 

where n is the density of the shocked gas, m e is the 
mass of the electron, o is the Steffan-Boltzmann constant, 
k is the Boltzman constant and e a number smaller than 1 
than accounts for the fact that the transfer of energy from 
the electrons to the dust grains is not completely efficient 
(eq. 4.27, Evans 1994). The shocked gas can reach tempera- 
tures as high as T 9as =10 6 K (Dopita & Sutherland 1996). 
Assuming n=10 4 ' 5 cm -2 (see §2.2.2), the temperature of the 
dust would be ^30 K. Since e must be <1, the heating we 
have calculated is an upper limit and therefore, the heating 
produced by collision with the electrons is not important. 

A more efficient heating source is the UV radiation from 
the shocked gas. Evidence for this radiation in NGC1068 has 
been discussed by several authors. Axon et al. (1998) found 
an increase in excitation of the gas (despite an increase in 
density) and an excess of continuum along the radio jet. The 
authors suggest that in addition to the central AGN, a local 
ionization source is present, which is likely to be the ionizing 
radiation from the shocked gas. 

We have checked whether this mechanism can explain 
the extended MIR emission in NGC1068. In order to do this, 
we have calculated the UV luminosity required to account 
for the MIR dust emission measured in the filters used by 
ALL00 and compared with the luminosity produced by the 
shocks (taking into account the higher absorption coefficient 
of dust grains at shorter wavelengths and the strong UV con- 
tinuum expected from incident radiation, we have assumed 
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Location 


t MIR 
^20 


t MIR 
^11 


T UV 

^20 


T UV 


T UV 

shock 


b (N-E) 
c (S-W) 


2.3 
2.1 


6.6 
0.8 


>57.5 
>52.5 


>253.0 
>30.7 


<106 
<130 



Table 3. Comparison between the UV luminosities required to 
sustain the observed MIR luminosities and the UV luminosity 
generated by the shocked hot gas. All luminosities are given in 
units of 10 41 erg s" 1 . None of the luminosities are monochromatic, 
but the luminosities in the filters discussed in the text. L|q and 
are the required UV luminosities to sustain the observed 
MIR luminosities Ljjj* and h^ IR , for the conversion efficiency 
calculated in the text. Lv > is the UV luminosity produced by 
the shocks (v s =700 km s -1 ). 



that the heating is dominated by photons with wavelength 
in the spectral range 30-3000 A , Ryter 1996). 

Only a fraction / of the incident radiation is absorbed 
by the dust. This radiation is reemitted in the infrared and 
only a fraction e is in the spectral window of a given filter. 
Therefore, the fraction of incident energy transformed into 
thermal energy in the filter is: 

V = e / 

Following the procedure by Alonso-Herrero et al. 
(1998), we obtain rj <0.04 (since / <1) for the filter with 
A o =20.5 /im and FWHM=1.1 /im, that is, 4% or less of 
the incident radiation is re-emitted as thermal dust emission 
within the 20.5 /im filter. We obtain rj <=0.026 for the filter 
with Ao=11.2 pm and FWHM=0.44 /im. For these calcula- 
tions we have assumed that the MIR emission in a given 
filter is dominated by dust with temperature such that the 
maximum flux is emitted at the central wavelength of the 
filter; i.e. 145 and 260 K respectively (ALL00 et al. conclude 
that the expected temperature is 150 in region c). 

The UV luminosities (Lij , L^) needed to sustain the 
MIR luminosities measured in regions b and c (Li{ IR , L2o IR ) 
are shown in Table 3 (Li{ IR , L2o IR have been calculated 
with the flux values given by AL00 and the properties of the 
two filters). Do shocks generate luminosities of this order? 

The UV ionizing radiation field (over 4-7T sr) produced 
by a shock of velocity v s propagating in a medium of (un- 
shocked) density n is given by eq. 4.1 in Dopita & Sutherland 
1996: 



rtZck = i.n x io" 3 



i3.04 n 



-2 -1 

erg cm s 



100 km s 1 cm 
For v a =700 km s" 1 and n =100 cm" 3 (see §2.2.2) the 
total UV luminosity emitted by the shocked gas in regions 
b and c is given in Table 3. To calculate the area occupied 
by the shocked gas, we have assumed that b and c have 
cylindrical geometry. The height and diameter of the cylin- 
der were taken as the projected sizes of regions b and c in- 
dicated in Table 1. We obtain Areab= 2.6xl0 41 cm" 2 and 
Area c ~ 3.2 x 10 41 cm" 2 (projected). The results of these cal- 
culations (Table 3) show that the shocks produce an intense 
continuum that is likely to contribute to the dust heating. 
According to these results, the shock continuum is enough 
to explain the MIR dust emission in c and the 20pm emis- 
sion in b. However, our calculations assume that the whole 
area of the cylinder is filled with shocked gas, which is un- 
likely and that all the photons are intercepted by the dust. 
This is certainly not the case, since part of the photons ion- 



ize the gas and part escape unabsorbed (Axon et al. 1998). 
Therefore, L^ ock is an upper limit and an additional heating 
source is probably present. This is actually needed to explain 
the 11/im flux measured in b. The active nucleus can pro- 
vide that extra heating. In order to subtain the 11/xm flux 
in region b we need L uv >7xl0 44 erg s _1 from the active 
nucleus. Using the information on the observed total spec- 
trum of NGC1068 (Pier et al. 1994), Marco & Alloin (2001) 
obtained a lower limit of L >4xl0 44 erg s _1 . Therefore, 
the AGN is likely to provide the additional heating source. 

We conclude that the shocked gas generates an intense 
UV continuum that enhances the MIR emission in the re- 
gions overlapping with the radio structures. The AGN con- 
tinuum also contributes to the dust heating. 

Another possibility mentioned in the previous section 
is that the radio jet has exposed the inner dusty cores of 
the clouds. In this case, the UV photons must be generated 
elsewhere, in the much faster moving emission line gas. The 
calculations done above are also valid in this case and both 
the shock and the AGN continuum are likely to contribute 
to the heating of the dust in the cores of the clouds. 



3 PKS2152-69 (Z = 0.028) 

The interaction between the radio jet and a high ionization 
cloud at a projected distance of 8 kpc from the nucleus has 
been widely studied (e.g. Tadhunter et al. 1987; Fosbury et 
al. 1998). The cloud lies approximately along the radio axis 
and is believed to be hit by the radio jet. The existence of 
dust in the cloud is confirmed by the detection of a very blue 
optical to near-UV continuum with 10 per cent polarization 
and the E- vector perpendicular to the position angle of the 
nucleus-high ionization cloud axis (di Serego Alighieri et al. 
1988). The authors explained the continuum and spectropo- 
larimetric properties of the cloud as the result of scattering 
by dust in the cloud of beamed radiation from the nucleus. 

Therefore, this appears as a clear example where shocks 
(produced by jet/cloud interactions) do not prevent dust 
from surviving in sufficiently large quantities, so it can still 
scatter the nuclear light and cause the alignment between 
the optical and radio structures. 

This does not imply, however, that shocks have not des- 
troyed the dust. High resolution radio maps and HST op- 
tical images (Fosbury et al. 1998) have shown that part of 
the high ionization cloud overlaps with a bright radio com- 
ponent. The optical continuum emitted by the overlapping 
region is redder and fainter than the continuum emitted by 
the rest of the cloud. This suggests that the dust grains 
responsible for the scattering have different properties. The 
change of continuum colours could be explained if the small 
dust grains (easier to destroy) have been destroyed by the 
shocks in the region of the cloud that is interacting with 
the radio structure. The scattering would be dominated by 
large grains, which produce a redder reflected continuum. 
Small grains dominate the scattering in the rest of the cloud 
(as proposed by Magris, Binette & Martin 1992) and there- 
fore the scattered continuum is bluer. Spatial spectropolari- 
metric information across the cloud and/or the study of the 
continuum colour variation would confirm whether this in- 
terpretation is valid or another explanation is required. 

In summary, PKS2152-69 is a radio galaxy that 
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presents close alignment between the optical and radio struc- 
tures. The aligned continuum is due to radiation from the 
active nucleus scattered by dust, in spite of the strong in- 
teraction between the radio structures and the cloud. In this 
case, the observations suggest that at least the small grains 
have been destroyed as a consequence of the interaction be- 
tween the gaseous cloud and the radio jet. The dominant 
source of aligned continuum is the region of the cloud that 
apparently has not been shocked. 

3.1 Discussion 

The case of NGC1068 implies that dust can survive under 
strong shock conditions, contrary to what we expect from 
our modeling. Possible explanations require recent interac- 
tions or dense molecular dusty cores exposed to the AGN 
radiation field. The close correspondence between the radio 
and dusty structures in this galaxy supports the model by 
Bremer, Fabian & Crawford (1997) for the alignment effect 
in high redshift radio galaxies as a consequence of scattered 
light. In this model the radio jet has exposed the inner dusty 
cores of the clouds and the contribution of scattered light is 
enhanced. An alternative possibility is that shocks heat the 
dust radiatively and enhance the MIR emission along the 
jet. 

The case of PKS2152-69 suggests that shocks can indeed 
destroy small dust grains (although spatially resolved spec- 
tropolarimetric observations of the high ionization cloud are 
needed to confirm this). However, dust survives in enough 
quantities in adjacent (or distant) regions to produce the 
apparent alignment between the radio and optical struc- 
tures. 

Both cases show that the interaction between the ra- 
dio jet and the ambient gas has not prevented dust from 
surviving in enough quantities to scatter light from the ac- 
tive nucleus and produce alignment between scattered light 
and the radio structures. This could also be the case in high 
redshift radio galaxies. 



4 SEARCHING FOR DUST IN SHOCKED 
REGIONS 

Another way to investigate the survival of dust in shocked 
regions is by studying the line emission from species sensitive 
to depletion, such as calcium and iron. In the ISM 99% of 
calcium and iron are locked in the dust grains, rather than 
in the gaseous phase. 

Villar-Martin et al. (1996, 1997) (based on the idea 
proposed by Ferland 1993) showed that in the ab- 
sence of dust, the doublet lines [CaII]A7291,7324 should 
be easily detectable (stronger than [OI]A6300) in pho- 
toionised nebula around nearby radio galaxies (such that 
[OIIIJA5007/H/3 <10). The weakness (or absence) of these 
lines implies the existence of dust mixed with the ionised 
gas if the dominant ionization mechanism is AGN photoioni- 
zation, whether shocks are present or not. This diagnostic 
method is not reliable when applied to AGN photoionised 
nebulae at higher redshift (z >1). There is evidence (Gond- 
halekar, Gondhalekar & Wilson 1988; Francis 1993) that the 
AGN continuum is harder at higher redshifts. With the ex- 
pected power law index (a—-l, Villar-Martin, Tadhunter & 



Clark 1997 instead of q=-1.5 as assumed for low redshift ra- 
dio galaxies) the [Call] emission lines should be very weak, 
unless for low ionization parameters. The non detection of 
[Call] emission at high redshift is therefore no so straight 
forward to interpret in terms of depletion by dust. 

Over the last few years, the dominant emission line 
mechanism (AGN photoionization vs. shock related pro- 
cesses) in high redshift radio galaxies has been a source of 
debate (e.g Villar-Martin et al. 1997, Best, Rottgering & 
Lehnert 2000; de Breuck et al. 2000). Both mechanisms are 
likely to play a role. 

We investigate here the applicability of the calcium 
test as a diagnostics for the existence of dust under shock 
ionization conditions. We have run a set of shock and 
shock+precursor models with the MAPPINGS III code, 
which is an enhanced version of the MAPPINGS II code 
(Dopita & Sutherland 1996). MAPPINGS III is designed 
to handle higher temperatures (up to 10 8 K) and hence 
faster shocks. The models considered here range from 250 
to 800 km s^" 1 which are appropriate for the physical con- 
ditions of gaseous nebula around radio galaxies (the maxi- 
mum velocity considered by Dopita & Sutherland was 500 
km s _1 ). We have predicted the strength of the [CaII]A7291 
emission line relative to other lines in the absence of dust. 
The abundances are solar and non depleted by dust. In Fi- 
gure 1 we plot the variation of the [CaII]A7291 flux (the 
strongest line in the doublet) relative to other emission lines 
with shock velocities (250-800 km s _1 ). The two sequences 
(solid and dashed lines) correspond to magnetic parameters 
and 4 respectively. The diagram on the left shows the 
pure shock models and the diagram on the right shows the 
shock+precursor models. The models show that over this ve- 
locity range [CaII]A7291 should be stronger that [OIII]A4363 
and HeIIA4686 and, for velocities below 500 km s _1 , the line 
should be stronger or comparable to [OI]A6300. 

I.e., the [CaII]X7291 line should be detectable in the 
extended gas of radiogalaxies (and the NLR of Seyferts) 
with shocks (whatever the ionization mechanism is), stronger 
than HeIIX4686 and [OIII]\4363 (and similar to or stronger 
than [OIJX6300). The absence or weakness of the line would 
imply that dust is mixed with the shocked gas and the shocks 
have not destroyed the dust. The applicability of this test to 
high redshift radio galaxies where AGN photoionization is 
the dominant ionization mechanism is not obvious, except 
for low values of the ionization parameter. 

This test was applied to the high ionization cloud in 
PKS2152-69 (Villar-Martm et al. 1997). The lack of [Call] 
emission was interpreted by the authors as evidence for the 
existence of dust mixed with the ionized gas. This adds ev- 
idence to the survival of dust in the cloud. 

An interesting case to compare the models with is the 
radio galaxy PKS2250-41 (z =0.31), where it has been pos- 
sible to isolate the shocked and the non shocked gas. Clark 
et al. (1997) found clear evidence for a strong interaction 
(and shocks) between the radio jet and the extended gas 
of this radio galaxy. Villar-Martin et al. (1999) made the 
spectral decomposition of the emission lines emitted by the 
interaction region. This led to the kinematic resolution of the 
shocked gas (FWHM 600-900 km/s) and the non perturbed 
gas (FWHM~150 km/s) in the main emission lines. The line 
ratios measured for the shocked gas were consistent with 
shock ionization models with velocities ~200-300 km s _1 . 
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Velocity Velocity 
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Figure 1. Variation with shock velocity (in km s" 1 of the flux of [CaII]A7291 relative to HeIIA4686, [OIII]A4363 and [OI]A6300 for 
shock models with no (left) and with (right) precursor contribution. Each line in the diagram represents the ratio Magnetic parameter 
~Yp2 ls 1 (dashed lines) and 4 (solid lines) muQ cm -3 / 2 . The shadowed region shows the velocity range of the models that reproduced 
the optical line ratios of PKS2250-41. For all models [CaII]A7291 should be stronger than [OIII]A4363 and HeIIA468. The absence of 
weakness of the line would imply that the shocked gas is mixed with dust and therefore, the shocks have not destroyed the dust. 



As Fig. 1 shows, if the shocked gas is dust free, [CaII]A7291 
should be stronger than HeIIA4686, [OIII]A4363 and similar 
or stronger than [OI]A6300. Since these two lines have been 
detected in the shocked gas of PKS2250-41 (Villar-Martm 
et al. 1999), broad [CaII]A7291 (FWHM~600-900 km ~ x ) 
should also be detected. If this is the case, we shall conclude 
that the gas is dust free (calcium is not depleted) and shocks 
have destroyed the dust. The non detection of the line will 
imply that the shocked gas contains dust. 



5 SUMMARY AND CONCLUSIONS 

The goal of this paper has been to study the survival of 
dust in low redshift active galaxies where shocks have been 
induced by the passage of the radio jet and link the con- 
clusions with the alignment effect in high redshift radio ga- 
laxies. 

The close correspondence between the extended MIR 
and radio structures in NGC1068 suggests that dust exists 
in the regions shocked by the radio jet. We show that, in 
addition to the AGN radiation, the continuum generated by 
the shocks contributes to the dust heating in these regions, 
enhancing in this way the MIR emission. 

Models appropriate for NGC1068 show that small dust 
grains should have been destroyed in less than 1000 years. 
We propose several explanations for the dust survival, such 
as a recent interaction. Another possibility is that the clouds 
consist of an external shell of ionised gas and dense mole- 
cular dusty cores. The radio jet has exposed the cores to 
the intense UV radiation field of the nucleus and the MIR 
emission is enhanced along the radio axis. This model was 
proposed by Bremer, Fabian & Crawford (1997) to explain 
the enhancement of scattered light along the radio axis in 
high redshift radio galaxies. 



The close radio/optical alignment in the radio galaxy 
PKS2152-69 (due to dust scattered light) and the evidence 
for a strong interaction between the radio jet and the gas, 
prove that jet induced shocks and radio/optical alignment 
due to dust scattered light are compatible. Based on the 
different continuum colour observed between the shocked 
and non shocked regions, we suggest that small dust grains 
have been destroyed in the region of the interaction. This 
has not prevented the alignment between the dust scattered 
continuum and the radio structures to be observed. 

In summary, there can be high redshift radio galaxies 
where the cloud properties are such that dust destruction 
is prevented (for instance dense dusty cores) or where the 
interaction happened very recently and there has not been 
enough time for the dust to be destroyed. Another possibility 
is that, as in PKS2152-69, dust has been destroyed in the 
shocked gas, but regions adjacent (or distant) to the radio 
structures contain enough dust to account for the observed 
alignment. 

Our conclusion is that jet/cloud interactions and the 
alignment effect due to dust scattered light are compatible. 

We propose to use the [CaII]A7291 line to investigate 
the existence of dust mixed with the shocked gas in active 
galaxies (radio galaxies, Seyferts) where shock ionization 
(with or without precursor contribution) or AGN photoioni- 
zation dominate the emission line processes. In the absence 
of dust, [CaII]A7291 from the shocked gas should be stronger 
than [OIII]A4363 and HeIIA4686 (both lines often detected). 
The absence of the [Call] line would suggest that the calcium 
is depleted onto dust grains in the shocked gas. PKS2250- 
41 appears as a good candidate to apply this test, because 
the shocked and the non perturbed gas have already been 
resolved kinematically. If shocks have destroyed the dust, 
broad (FWHM~600-900 km s" 1 ) [CaII]A7291 should be de- 
tected. 



© 2001 RAS, MNRAS 000, 



8 M. Villar- Martin et al. 



6 ACKNOWLEDGMENTS 

We thank an anonymous referee for useful comments that 
helped to improve the paper. We thank A. Wilson for pro- 
viding the radio map of NGC1068. MVM thanks Carlos de 
Breuck for useful discussions. The work of LB was supported 
by the CONACyT grant 32139-E. 



REFERENCES 

Alloin D., Pantin E., Lagage P., Granato C, 2000, A&A, 363, 926 
(ALL00) 

Alonso-Herrero A., Simpson C, Ward M., Wilson A., 1998, ApJ, 
495, 196 

Axon D., Marconi A., Capetti A., Macchetto F., Schreier E., 

Robinson A., 1998, ApJ, 496, L75 
Best P., Rottgering H., Lehnert M., 2000, MNRAS, 315, 21 
Bicknell G., Dopita M., O'Dea C., 1997, ApJ, 485, 112 
Bremer M., Fabian A., Crawford C, 1997, MNRAS, 284, 213 
Capetti A., Axon, D., Macchetto F., Marconi, A., Winge, C, 

1999, ApJ, 516, 187 

Chambers, K. C, Miley, G. K., & van Breugel, W. 1987, Nature, 
329, 604 

Cimatti A., Dey A., van Breugel W., Antonucci R. , Spinrad R. , 

1996, ApJ, 465, 145 
Clark N., Tadhunter C.N., Morganti R., Killeen N.E.B., Fosbury 

R.A.E., Hook R.N., Siebert J., Shaw M.A., 1997, MNRAS, 

286, 558 

Draine B., Lee H., 1984, ApJ, 285, 89 

De Breuck C, Rttgering H., Miley G., van Breugel W., Best P., 

2000, A&A, 362, 519 

De Young D.S., 1998, ApJ, 507, 161 

Di Serego Alighieri S., Courvouisier T., Fosbury R., Tadhunter 

C, Binette L., 1988, Nature, 334, 591 
Dopita M. & Sutherland R., 1996, ApJS, 102, 161 
Evans A., in 'The dusty Universe', Ellis Horwood Series in Space 

Science and TechnologyPraxis Publishing LTD. 
Ferland G., 1993, in Proc. Madrid Meeting on The Nearest Active 

Galaxies, ed. J.E. Beckman, H. Netzer & L. Colina, p. 75 
Forman W., Jones C, Tucker W., 1985, ApJ, 293, 102 
Fosbury R., Morganti R. , Wilson W., Ekers R., di Serego Alighieri 

S., Tadhunter C, 1998, MNRAS, 296, 701 
Francis P., 1993, ApJ, 407, 519 

Kellermann K., Pauliny-Toth I., Williams P., 1969, ApJ, 157, 1 
Leep S., McCray R., Shull J., Woods D., Kallman T., 1985, 288, 
58 

Marco O., Alloin D., 2001, A&A, 353, 465 
Magris G., Binette L., Martin P., 1993, Ap&SS, 205, 141 
McCarthy, P. J., Spinrad, H., Djorgovski, S., et al. 1987, ApJ, 
319, L39 

Netzer H., Turner T., 1997, ApJ, 488, 694 

O'Brien P., Gondhalekar P., Wilson R., 1998, MNRAS, 233, 601 
Packham C, Young S., Hough J., Axon D., Bailey J., 1997, MN- 
RAS, 288, 375 

Pier E., Antonucci R., Hurt T., Kriss G., Krolik J., 1994, ApJ, 

428, 124 
Rytcr 1996, ApSS, 236, 285 
Scheuer P., 1995, MNRAS, 277, 331 
Seyfert C, 1943, ApJ, 97, 28 

Shirey R., Soria R., Borozdin K., et al., 2000, A&A, 365, L195 
Spitzer L., 1941, ApJ, 93, 369 

Stocke J., Rieke G., Lebofsky M., 1981, Nature, 294, 319 
Tadhunter C, Fosbury R., Binette L., Danzinger I., Robinson A., 

1987, Nature, 325, 504 
Tadhunter C, Scarrot S., Rolph C, 1990, MNRAS, 246, 163 
Tadhunter C, 1991, MNRAS, 251, 46 
Ulrich et al. 1981, MNRAS, 192, 561 



Vernet J., Fosbury R., Villar-Martm M., Cohen M., Cimatti A., 

di Serego Alighieri S., Goodrich R., 2001, A&A, 366, 7 
Villar-Martm M., Binette L., 1996, A&A, 309, 97 
Villar-Martm M., Binette L., 1997, A&A, 317, 350 
Villar-Martm M., Tadhunter C, Clark N., 1997, 323, 21 
Villar-Martm M., Tadhunter C, Morganti R., Clark N., Killeen 

N., Axon D., 1998, A&A, 332, 479 
Villar-Martm M., Tadhunter C, Morganti R., Axon D., Koeke- 

moer A., 1999, MNRAS, 307, 24 
Westerlund B. & Smith L., 1966, Australian J. Phys, 19, 181 
Wilson A. & Ulvestad J., 1987, ApJ, 319, 105 



© 2001 RAS, MNRAS 000, 0-?? 



